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We have developed an exact, general method to compute Casimir interactions between a finite 
number of compact objects of arbitrary shape and separation. Here, we present details of the 
method for a scalar field to illustrate our approach in its most simple form; the generalization to 
electromagnetic fields is outlined in Ref. [1]. The interaction between the objects is attributed to 
quantum fluctuations of source distributions on their surfaces, which we decompose in terms of mul- 
tipoles. A functional integral over the effective action of multipoles gives the resulting interaction. 
Each object's shape and boundary conditions enter the effective action only through its scattering 
matrix. Their relative positions enter through universal translation matrices that depend only on 
field type and spatial dimension. The distinction of our method from the pairwise summation of 
two-body potentials is elucidated in terms of the scattering processes between three objects. To 
illustrate the power of the technique, we consider Robin boundary conditions (j) — ^dn4> = 0, which 
interpolate between Dirichlct and Neumann cases as A is varied. We obtain the interaction between 
two such spheres analytically in a large separation expansion, and numerically for all separations. 
The cases of unequal radii and unequal A are studied. We find sign changes in the force as a function 
of separation in certain ranges of A and see deviations from the proximity force approximation even 
at short separations, most notably for Neumann boundary conditions. 

I. INTRODUCTION 

Casimir forces arise when the quantum fluctuations of a scalar, vector, or even fermion field are modified by the 
presence of static or slowly changing external objects [2] . The objects can be modeled by boundary conditions that they 
place on the fluctuating field by an external fleld, cr, to which (j) couples [3], or, in the case of electromagnetism, 
by a material with space and frequency dependent dielectric and magnetic properties. The Casimir energy is the 
difference between the energy of the fluctuating field when then objects are present and when the objects are removed 
to infinite separation. 

The advent of precision experimental measurements of Casimir forces [4] and the possibility that they can be applied 
to nanoscalc electromechanical devices [5, 6] has stimulated interest in developing a practical way to calculate the 
dependence of Casimir energies on the shapes of the objects. Many geometries have been analyzed over the years, but 
the case of compact objects has proved rather difficult. In a recent Letter [1] we described a new method that makes 
possible accurate and efficient calculations of Casimir forces and torques between any number of compact objects. 
The method applies to electromagnetic fields and dielectrics as well as perfect conductors. It also applies to other 
fields, such as scalar and Dirac, and to any boundary conditions. In this approach, the Casimir energy is given in 
terms of the fiuctuating field's scattering amplitudes from the individual objects, which encode the effects of the shape 
and boundary conditions. The scattering amplitudes are known analytically in some cases and numerically in others. 
If the scattering amplitudes are known, then the method can be applied from asymptotically large separation down 
to separations that are a small fraction of the dimension of the objects. Results at large separations are obtained 
using low frequency and low angular momentum expansions of scattering amplitudes. The coefficients multiplying 
the successive orders in inverse separation can be identified with increasingly detailed characteristics of the objects. 
At small separations the manipulation of large matrices, whose dimensions grow with angular momentum, eventually 
slows down the calculation. However at these distances other methods, notably the "proximity force approximation" 
(PFA) , apply. Thus it is now possible to obtain an understanding of Casimir forces and torques at all separations for 
compact objects. 

The aim of this paper is to provide a pedagogical introduction to our methods by treating in detail the simplest case, 
a scalar field obeying a boundary condition on a sharp surface. The complications of electromagnetism and smoothly 
varying dielectrics were already introduced briefly in Ref. [1]. They will be treated in more detail in subsequent 
publications [7]. Our approach relies on a marriage of methods from path integral and scattering formalisms, so we 



2 



provide background on both of these subjects as they apply to Casimir effects. 

Schwinger, in particular, emphasized that Casimir forces could be understood as ordinary electromagnetic interac- 
tions between quantum fluctuations of charge and current in metals or dielectrics [8]. We implement this idea here 
through a functional integral formulation. We begin by writing the Casimir energy as the logarithm of a functional 
integral over all field fiuctuations constrained by the boundary conditions on a set of surfaces. Following Refs. [9, 10] 
we implement the boundary conditions by introducing integrals over sources that enforce the boundary conditions 
as functional (5-functions. Next we perform the functional integral over the field. The result is an integral of the 
form Ha/ T^Qaexp{iS[Q]), where Qa are the sources on the different surfaces labeled by a and S[g] is the classical 
action of the field, which is uniquely determined by the sources g^. The contributions to S[g\ are of two qualitatively 
different forms, those involving sources on different surfaces and those coupling a particular surface source to itself 
through the field it generates. Both have simple representations in a basis of angular momentum eigenstates ("partial 
waves") and multipole moments of the sources. The couplings between different objects involve only a well-known 
kinematic "translation matrix," U, that relates a partial wave amplitude generated by one source to partial waves seen 
by another [11]. The self-interactions depend only on the "transition matrix," T, (related to the scattering matrix, S, 
by T = i(S — I)]) [12], which describes the response of the scalar field to the boundary condition on the surface. The 
result is remarkably simple. For a complex scalar field in the presence of two objects it takes the form. 



where the determinant is over the partial wave indices on the matrices T" and U"'^ and the integral is over k = —iui/c, 
the imaginary wavenumber. In Section IV the usefulness of this result is demonstrated through several specific 
applications. 

Casimir forces between compact objects were first considered by Casimir and Polder in 1948 [13]. Since then, two 
threads of work related to ours have been pursued: attempts to evaluate the Casimir force between compact objects 
explicitly, and efforts to develop a general framework similar to that embodied in Eq. (I.l). Until recently work along 
the first line consisted of expansions at asymptotically large separation. Casimir and Polder found the electromagnetic 
force between two neutral but polarizable atoms to leading order at large separation [13]. Feinberg and Sucher [14] 
generalized to arbitrary compact objects and included magnetic effects. Balian and Duplantier studied perfect metals, 
and derived explicit results to leading order at asymptotically large separation [15]. More recently Gies et al. [16] 
used numerical methods to evaluate the Casimir force between two Dirichlet spheres for a scalar field, over a range 
of subasymptotic separations, and in other open geometries such as a plate and a cylinder [17] or finite plates with 
edges [18]. Bulgac and collaborators [19] applied scattering theory methods to the same scalar Dirichlet problem and 
obtained results over a wide range of separations. The only explicit calculations for subasymptotic distances up to 
now have been for a scalar field obeying Dirichlet boundary conditions on two spheres, a sphere and a plate [19] and 
for electromagnetic fields for a plate and a cylinder [20] and two perfectly conducting spheres [1] . 

Formulas for the Casimir energy closely related to Eq. (I.l) have appeared previously in the hterature. The first 
appearance we are aware of is in the work of Balian and Duplantier [15] based on the multiple refiection expansion 
(MRE). Their Eq. (7.20) gives a quantity, that is directly related to the Casimir energy, and could form the basis 
for an approach like ours. Their approach, like ours, punctuates periods of propagation in the vicinity of each object 
with free propagation between the objects. However no explicit expressions for the propagation kernels, Ka, are given. 
More recently "log-det" formulas like ours have arisen in several functional integral studies of Casimir forces [20, 21]. 
To our knowledge Kenneth and Klich in Ref. [22] were the first to identify the inverted Green's function as a T-matrix 
and derive Eq. (I.l) as a formal result. Here also, no explicit form for the matrices was derived. Their derivation 
applies to scalar fields in a medium with a space and frequency dependent speed of light. For the case of a scalar 
field, our result can be viewed as an explicit expression for their formula in a basis of partial waves, which we show 
is particularly well suited to practical applic;ations. Our new derivation of Eq. (I.l) in terms of quantum fiuctuations 
of sources has the advantage that it allows for a physically transparent extension to gauge fields in the presence of 
material objects with general dielectric and magnetic properties, as outlined in Ref. [1]. 

In Section IV, we provide a number of examples where our approach yields straightforward results for Casimir 
interactions that could not be obtained before. We carry out explicit analytical and numerical computations for 
spheres with Robin boundary conditions (f) — Xdn(t> = 0. The sign of the Casimir interaction at asymptotically 
large and small distances is classified as a function of the value of the Robin parameter A. The Casimir energy at 
asymptotically large separations is obtained as a series in the ratio of distance to sphere radius. We also perform 
numerical computations of the Casimir interaction of two spheres with Robin boundary conditions spanning the 
full range of separations, including very short distances, where our results confirm the proximity force approximation 
(PFA). At intermediate distances considerable deviations from the PFA are found. For certain values of A, the Casimir 
force can change sign (once or twice) as function of separation. 




(1.1) 
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Finally, in Sec. IV D, wc give an expression for the Casimir interaction at large separations between two compact 
object of arbitrary shape and with arbitrary boundary conditions in terms of generalized capacitance coefficients. 
Technical derivations are presented in three Appendices. 



II. FOUNDATIONS 



In this Section we review formalism essential for our work. We use both functional integral methods and techniques 
originating in scattering theory. First, in Sec. II A we introduce the functional integral approach of Refs. [10, 23, 24], 
which allows us to trade the problem of fields fluctuating in the bulk for the interactions of sources defined only on 
the bounding surfaces. In our approach, individual objects are characterized by the way they scatter the fluctuating 
fields. This information is summarized in the transition matrix. T. which we introduce in Sec. II B and relate to 
scattering solutions of the equations of motion and to Green's functions. In the A; ^ limit, Helmholtz's equation 
reduces to Laplace's equation. In this limit the T-matrix can be related to generalized coeflicients of capacitance, 
which we summarize in Sec. II C. Finally in Sec. II D we introduce the translation formulas that relate partial waves 
computed with respect to one origin to those computed with respect to another. 



A. Functional integral formulation 

We consider a complex quantum field, 4>{x, t), which is defined over all space and constrained by boundary conditions 
C on a set of fixed surfaces E^, for a = 1, 2, . . . , iV, but is otherwise non-interacting. We assume that the surfaces are 
closed and compact and refer to their interiors as "objects." Our starting point is the functional integral representation 
for the trace of the propagator, Tr e"*-^'^-^/^ [25], 



j [D<j)\c eH^^I^I = Z[C], (II. 1) 



where the subscript C denotes the constraints imposed by the boundary conditions.^ The integral is over all field 
configurations that obey the boundary conditions and are periodic in a time interval T. S[(t)] is the action for a free 
complex field. 



Se[<I>] = - l^-^l") ' ("-2) 

where the x-integration covers all space. ^ 

The ground state energy can be projected out of the trace in Eq. (II. 1) by setting T = —iA/c taking the limit 
A — > oo, 

£o[C] = - lim ^ In f Tr e"^^^/''^') = - lim ^ In Z\C] , (II.3) 

and the Casimir energy is obtained by subtracting the ground state energy when the objects have been removed to 
infinite separation, 

he 

£:[C] = - lim ^ln(Z[C]/Z^) . (II.4) 

A— >cx) A 

In the standard formulation, the constraints are implemented by boundary conditions on the field ^ at the surfaces 
{Eq}. The usual choices are Dirichlet, (/) = 0, Neumann, 9„(/) — 0, or mixed (Robin), (\)—\dn<\) = 0, where 9„ is the 
normal derivative pointing out of the objects. To be specific, we first consider Dirichlet boundary conditions. The 
extension to the Neumann case is presented in Appendix A. As noted in the Introduction, the only effect of the choice 
of boundary conditions is to determine which T-matrix appears in the functional determinant, Eq. (I.l). 



^ We have used an abbreviated notation for the functional integral. Since </> is complex / T>if> should be understood as / T>if>T>if>* , and 
similarly in subsequent functional integrals. 

^ Note that <j> is defined and can fluctuate inside the objects bounded by the surfaces S^. In this feature our formalism departs from 
some treatments where the field is defined to be strictly zero (for Dirichlet boundary conditions) inside the objects. The fluctuations 
interior to the objects do not depend on the separations between them and therefore do not affect Casimir forces or torques. 
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For readers who are not familiar with the functional integral representation of the Casimir energy, Eq. (11.4), in 
Appendix B we show that Eq. (II. 4) agrees with the traditional definition of the Casimir energy in terms of zero-point 
energies of normal modes. 

Since the constraints on (/> are time independent, the integral over (j){x, t) may be written as an infinite product of 
integrals over Fourier components, 



J mc = n \.'^M^)]c ' 



(II.5) 



where 



0(x,i)= J2 <Pn(x)e^™*/^, 



(II.6) 



and the logarithm of Z becomes a sum. 



lnZ[C]= f2 ln|y[I?<^n(x)]cexp /'^'^ H ^ ) - l^M^Y^ 



(II.7) 



As T — ^ oo, J2n replaced by |^ dk, where k = 2im/{cT) and 0„(x) is replaced by 0(x, k). Combining 

the positive and negative fc-integrals gives 

hiZ[C] = (ifcln|y"[r»</)(x,fc)]cexp i'^ J dx{k^\(P{x,k)\^ - \V<p{x,k)\'^) | 

cT f°° 

= — / dk\n^c{k), (II.8) 



where 



3c(fc) = j [C</)(x,fc)]cexp j rfx(fcV(x,fc)|'-|V<A(x,A;)|2) 



(11.9) 



is the functional integral at fixed k. 

To extract the Casimir energy, we use T = —iA/c and Wick rotate the fc-integration (fc = m with k > 0). ^ Using 
Eq. (II.4), we obtain. 



Jo 3<x>(«k) 



Here 5c(*'«) is given by the Euclidean functional integral 

T 

h 



3c(«k) = j P'/'(x, iK)]^ exp 



dx(K^|0(x,m)p + |V0(x,iK)p) 



(11.10) 



(11.11) 



It remains to incorporate the constraints directly into the functional integral using the methods of Refs. [9, 10]. 
Working in Minkowski space, we consider the fixed frequency functional integral, 3c{k) (and suppress the label k on 
the field (p). Following Ref. [9, 10], we implement the constraints in the functional integral by means of a functional 
(5-function. For Dirichlet boundary conditions the constraint reads. 



j [I?<^(x)]c = j [Dct>{x)]l[ j [2^e«(x)]exp i^j^ dx(e*(x)<^(x) + c.c.) 



(11.12) 



A more careful treatment of the rotation of the integration contour to the imaginary axis is necessary in the presence of bound states. 
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where the functional integration over (f) is no longer constrained. Other boundary conditions can be implemented 
similarly. In the resulting functional integral, 



3c(fc) = n/ [DQai^)] I [m^)]e^p dK{k'\ct>i^)f-\V<f>{^)\')+J2l^ dx(^*(x).^(x)+c.c.) 



a=l 



n / [Ve^{^)] f [D.^(x)]exp 



T ~ 



(11.13) 



the fields fluctuate without constraint throughout space and the sources fluctuate on the surfaces. We denote 
the new "effective action" including both the fields and sources by S\(^, £»] . 



B. The scattering amplitude 

In this section we consider each of the objects in isolation, and review the solution to the Helmholtz equation 

- (V2 + 0^(x) = 0, (11.14) 

in the domain outside S^. We assume that (f)a obeys the Dirichlet boundary condition on S^. It is convenient to 

fix origin, C^a-) 

some point inside the object and introduce a coordinate vector, Xq,, defined with respect to this 
origin. A specific choice of the origin will be made later to simplify the analysis. For simplicity we suppress the label 
a from now on in this section. We introduce spherical polar coordinates relative to this origin, r = |x| and x. Outside 
E, ^ is a superposition of partial waves with definite angular momentum, 

</)(x) = J2 [cimhl'^ikr) + dimh'i'Jlikr)] YU^) , (11.15) 

Im 

where the spherical Hankel functions, h^i,^^{kr) and hl^,l{kr) are the outgoing (asymptotic to e^^^ jkr) and incoming 
(asymptotic to e~'^''^/kr) solutions respectively. To make contact with traditional methods of scattering theory, we 
rewrite ^{x) in terms of solutions with unit incoming amplitude, 

^i^) = m '^lrn<l>lm{^), (11.16) 
Im 

where 

0j™(x) = h'l^\kr)Yim{±) + Y^Si,m'im{k)h[}\kr)Yi,m'{^) ■ (11.17) 

The coefficients Simi'm'{k) measure the response of the object, S, to a unit amplitude incoming wave with angular 
momentum (Im). They are the matrix elements of the scattering operator, or §-matrix, 

Svm'imik) = {l'm'\S{k)\lm) , (11.18) 

which are fixed by the condition that the resulting basis functions, (j)im vanish on S. With the object absent, § would 
go to the identity operator, Si'm'im SwSmm', and (j)im{x) would reduce to 2ji{kr)Yim{i), the partial wave solution 
regular at the origin. It is convenient to make this result explicit by rewriting Eq. (11.17) as 

<^;„(x) = 2ji{kr)YU±) + Y,'^Ti,^,im{k)h^v\kr)Yi,^,{yi) , (11.19) 

I'm' 

where T is the transition matrix or T-matrix^ 

Ti,m'lm{k) = {l'm'\T{k)\lm) = ^ {Si^m'lmik) - Si^iSm'm) ■ (11.20) 



* Our definition of tfie T-matrix differs by a factor i from some conventional cfioices [26]. 
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Unitarity requires S'''S = I, or X^i'm' ^i'm'i"m"^i'ni'im = Si"i5m"m- Time reversal symmetry requires that Si'mUm is 
symmetric in {l,m) ^ {I' ,m'). If the object E is spherically symmetric, § and T are diagonal in £ and m. However 
we are also interested in more general cases where the material objects do not have any special symmetry. 

It is useful to review the connection between Ti>m'im{k) and the scattering of (j) from the object. The scattering 
amphtude is defined as the response of the object to an incoming plane wave. Since angular momentum is not 
conserved, we have to keep track of the direction of the incoming plane wave. The solution to the Helmholtz equation 
that reduces at r ^ oo to a plane wave with wave vector k accompanied by an outgoing scattered wave defines the 
scattering amplitude, /(k, k'), 

$(k,x) = 27r^0,„(x)z'y4(k) 

Im 

~ e''""+^/(k,k') asr^oo, (11.21) 
where k' = fc x is the wave vector of the observed wave. The differential cross section is given by 

| = i/(k,k')r (11.22) 

To relate /(k, k') to T we express the plane wave a a superposition of partial waves, 

e'k-x = 4nJ2Mkr)i'YUm*M , (11.23) 

Im 

where the unit vectors x and k specify the polar angles of x and k with respect to an arbitrary axis. We identify 
what remains with the scattered wave, 

$(k,x) = e'^-^ + An J2 i%m'im{k)h\}\kr)Yi,m'{mi*M- (n.24) 

I'm'lm 

Using the asymptotic form of the Hankel functions, h\^\kr) — ^ i~^~^e^^^ /kr^ we obtain 

/(k,k') = J E i'-''Tvm'imik)Yi,^'{^')YrM- (11.25) 

I'm'lm 

When angular momentum is conserved, Si'm'im{k) = e^^^^^^^^ii'^mm' ^tnd we obtain the standard result /(k, k') = 
Ej(2;+ l)/j(fc)P((cos0kk'), with fi{k) = isin<5,(fc)e^*'W. 
Finally, we will need the partial wave expansion for the free, outgoing wave Helmholtz Green's function, 

eo(x,x',/c) = ^TS^n = ^kY,3i{kr^)h^^\kr>)Yim{ii)Yil,{^) = ikJ2ji{kr<)h['\kr>)YUi')Yil{ic) , (11.26) 

Im Im 

where the notations r<(>) refer to whichever of r,r' is the smaller (larger). 



C. Low energy scattering and the coefficients of capacitance 

As fc — > 0, the Helmholtz equation reduces to Laplace's equation. Therefore we can relate the low energy limit of the 

§-matrix elements to the parameters that describe solutions to Laplace's equation, which are tensor generalizations of 
capacitance. In this subsection we continue to suppress the label a that distinguishes the particular object of interest. 

At large distances, solutions to Laplace's equation are of the form r'y;m(x) or r~'~^yim(x). If an object is placed 
in an external field (potential) of the first form, then its response, determined by the boundary condition (j) = Q on the 
surface, is of the second form. The monopole response to an asymptotically constant field defines the capacitance; the 
dipolc response to an asymptotically dipole field defines the polarizability, and so forth. Therefore we parameterize a 
solution to Laplace's equation that goes to r^Yim{^ by^ 

<^;„(X) OC r'yi„(x) - Cl'm'lm^^'^^ • (11.27) 
I'm' 



^ Note the minus sign that preserves the usual definition of capacitanee. When an object is held at a voltage V, its charge is Q = CV. 
Equivalently, when the object is grounded and the potential at infinity is V, the induced monopole field is —CV jr. 
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For a Dirichlct boundary condition the {Ciim'irn} are the tensor generahzations of the capacitance. For other boundary 
conditions, the physical connection to electrostatics is lost, but the tensor structure (and the connection to the 5- 
matrix) remains. 

Taking the limit k in Eq. (11.19) we obtain an expression for in terms of the k limit of the T-matrix, 

2(krV tr-^ (21' — 1 V 

I'm' 

In order to obtain a finite and non- vanishing limit, we find Ti/m'lm{k) ^ k^'^'' as fc — > 0, a standard result. Comparing 
these two expressions we find 

= (2/+"l)^?iS'-l)!! • ("-'^^ 

Some special cases deserve mention. For ? = /' = 0, as fc ^ 0, '7oooo(fc) ~ —ika, where a is the scattering 
length. Comparing Eqs. (11.28) and (11.29) we see that the capacitance, C = CqooO) equals the scattering length. The 
I — 0,1' = 1 coefficient describes the dipole moment developed by a conductor in response to a constant external 
potential (or, by symmetry, the charge induced by an external dipole field). It vanishes for spherically symmetric 
objects, but more generally we can choose the origin of the coordinate system so that Cimoo — Cooim = 0. From 
now on we assume that the origin, Oa, within the object S„ has been chosen to eliminate this dipole response. The 
I = I' = 1 coefficients describe the dipole response to an external dipole field: these nine components form the tensor 
(electrostatic) polarizability. 

Finally, we note that the coefficients Ci'm'im can be interpreted as the components of an abstract tensor C^'^ in a 
spherical basis. Cfi is the exterior product of two irreducible tensors, one rank I and the other rank I' . It is sometimes 
useful to project Cj'; into its irreducible components. Thus, for example, Cn, the electrostatic polarizability, can be 
decomposed into its trace, the scalar polarizability, and a traceless part that describes the aspherical response of the 
object. 



D. Translation formulas 



Partial wave solutions to the Helmholtz equation for an object Sq,, defined with respect to the origin O,,,, can be 
expanded in partial waves with respect to a second origin Op within a second object S^. The objects, coordinate 
origins, and notation for vectors are shown in Fig. 1. Any solution with definite angular momentum defined with 
respect to Oa, whether regular (ji) or outgoing {h'l^''), is regular when viewed from O/3. Therefore it is possible to 
expand both in terms of spherical Bessel functions jiikr^) with respect to Op, 

jiikra)YU^a) = E^;'m';m(X/3a)i/'(fcr-;3)y,,„,(x^) 

I'm' 

h'j'\kra)Yir^{±a) = Y.^it'imMji'ikrp)Y,^'iifi) (11.30) 

I'm' 

for rff < dap- Note that Xq, and Tip refer to the coordinate components of a point x relative to the origin Oa 
or Ojj respectively. The matrices U and V, that define the change of basis depend on the vector from Oa to Op, 

defined in Fig. 1. These formulas are known in the literature as translation formulas. 
A brief derivation of the translation formulas can be found in Appendix C and further discussion in Ref . [27] . 

Without loss of generality, to simplify the formulas we take the Cartesian unit vectors that define the orientations 
of the reference frames all to be parallel. We will need only formulas that apply to the outgoing solutions. 



Kt'lmi^cp) = VA^{-ir'i'-''^{2l + l){2l' + l)^^ 

X 5^ (-l)™"^'"^/2F^ (^J I ^ J^, ^'^,)jh\^){kdap)Yi„^"{±a0) (11.31) 

The summation over I" involves only a finite number of terms since the S-j symbols vanish for I" > l+l' and I" < \l^l'\. 
Moreover they are zero if I + I' + I" is odd or if ni" ^ m — m! . The matrix U'^^^n^ has the following symmetries, 

'^I'm'lm — \~^> '^Iml'm' ' 

^Vm'lm = ^V-m'l-m (11.33) 
^i^m'lm = i-^r'l^fm'lm- (11.34) 



^1 \Xi / 
\ / 

\ / 

V 

FIG. 1: (Color online) Two objects enclosed by the surfaces Ei and E2 are shown, each with bounding spheres (radii Ri 
and R2)- We assume that it is possible to choose bounding spheres that do not overlap. Coordinate systems with parallel 
axes are erected at origins Oi and O2 at positions Xi and X2. The origins are chosen so that the dipole coefficients of 
capacitance, Cooim, vanish. Coordinate vectors Xi and X2 to an arbitrary point, x, are shown. The vector from Oi to O2 is 
X12 — X2 — Xi = xi — X2. The distance between the two objects is di2 — |Xi2|. 

If we are considering only two objects, then it is always possible to orient the Cartesian coordinate systems so that 
they are related by translation along the z-axis. In this case simplifies. For Xq,^ — ±daf3Z one obtains, using 

Y,„„,„{±i) {±y"^{2U'+l)/4TT6,n"0, 

l^^L'imi±dc.pi) = <5™-™(-l)™*'- V(2^ + 1)(2Z' + 1) ^ i^''{2r + 1) (^J I o)4V(fc<ic.0). (11.35) 

III. EVALUATION OF THE CASIMIR ENERGY 

A. Performing the integral over 

We start with the expression for the fixed-A: functional integral, Eq. (11.13). For any fixed sources, {ga}, there 
is a unique classical field, ^iciip]! that is the solution to 6S[(p, g]/S(p{x.) = 0. The classical theory defined by S[(p, g], 
describes a complex scalar field coupled to a set of sources on the surfaces, and is a generalization of electrostatics. 
By analogy with electrostatics, the field is continuous throughout space, but its normal derivative jumps by ga{^) 
across S^. Indeed, the classical equations of motion that follow from SS/5(j) — are 

(V^ + ^^j(x) = 0, forx^S,, 
A(?!)ci(x) = 0, for X e Ea, 
AdnM^ = ffaW, forxeE^, (III.l) 

where A(p — (^in — (/)out and A9„(^ — dn4'\in — (?„0out- The subscripts "in" and "out" refer to the field inside and 
outside the bounding surface Sq,. As before, all normals point out of the compact surfaces. The solution to Eq. (III.l) 
is unique up to solutions of the homogeneous equations, which we exclude by demanding that (jjci vanish when the 
{ga} = 0. Continuing the analogy with electrostatics, we can write the classical field in terms of the free Green's 
function and the sources. 



'/'ci(x) = ^ / dx'5o(x,x',fc)£i;3(x') , 

(3 •''^13 



{111.2) 
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where Qo is the free Green's function given in Eq. (11.26). 

To compute the functional integral over </), we first decompose 
fluctuating part, 



into the classical part given by Eq. (III. 2) and a 



0(x) =,^cl(x)+(50(x). 

Then, because the effective action, S, is quadratic in <p, the 64> dependent terms are independent of (pd, 

N 



/ (ix(fc2|5,^(x)|2- |V(5<^(x)l2^ 



(III.3) 



(III.4) 



The classical action can be simplified by using the equations of motion, Eq. (III.l), which make it possible to express 
the action entirely in terms of integrals over the surfaces {Sa}) 



SciIq] = \Y.I (^a W-^cKx) + C.C. ) 



(III.5) 



where <?!)ci(x) is understood to be a functional of the sources 

The functional integral over 5(j) is independent of the classical field <j)c\ and defines the energy of the unconstrained 
vacuum fluctuations of (j). This term is divergent, or, more precisely, depends on some unspecified ultraviolet cutoff. 
However it can be discarded because it is independent of the sources and therefore common to 3c and 3oo- Note that 
this result is an explicit demonstration of the contention of Ref. [28]: the Casimir force has nothing to do with the 
vacuum fiuctuations of (j), but is instead a consequence of the interaction between fiuctuating sources in the materials. 
It is therefore not directly relevant to the fiuctuations that are conjectured to be associated with the dark energy. 

Prom Eq. (III. 2) it is clear that the solution to Eq. (III.l) obeys the superposition principle: ^!'ci(x) is a sum of 
contributions from each of the sources. 



0cl(x) 



/3 



(III.6) 



where (l)fj satisfies Eq. (III.l) with all sources set equal to zero except for gp. So the action can be expressed as a 
double sum over surfaces and over contributions to 4>c\ generated by different objects. This leaves a partition function, 
3c (^) 5 of the form 



N 

3c(fc)= n / [2^^a(x)]exp l\Y.j rfx(^;(x)(/,0(x) + c.c.) 

Q=l a, (3 



(III.7) 



to be evaluated. 



B. Evaluation of the Classical Action 



The classical action in Eq. (III. 7) contains two qualitatively different terms, the interaction between different sources, 
a ^ /3, and the self-interaction of the source Both can be expressed as functions of the multipole moments of the 
sources on the surfaces. 



1. Interaction terms: a ^ (3 



Consider the contribution to the action from the field, ^/j, generated by the source, q^^ integrated over the surface 

S0a = ^ j dx^ (£>* (Xc,)0/3(Xq,) + C.C. ) , (III.8) 

where the subscript on Xq indicates that the integration runs over coordinates measured relative to the origin of 

object a. The field (/)^(x^), measured relative to the origin of object /?, can be represented as an integral over its 
sources on the surface S/3 as in Eq. (III. 2). Since every point on is outside of a sphere enclosing S^, the partial 
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wave representation of Qo simplifies. The coordinate x> is always associated with X/3 and x< is identified with x^, so 
Eq. (III. 2) can be written 

M^0) = ikYl {kr0)YU±0) [ dx'^iKfcr-;3)F4(i'^)^/3(x'^) . (III.9) 

Im ''^P 

Note that the arguments of the Bessel functions and spherical harmonics are all defined relative to the origin Op. In 
particular, r'fj and xjlj are the radial and angular coordinates relative to Of^ corresponding to a point x' on the surface 
The integrals over define the multipole moments of the source g^, which will be our final quantum variables, 

Q0,im = / dx0ji{kr0)Yi*^{x0)Q0{-Kfs), (III.IO) 

so that 

'P/si'^fs) = ik^Qp,imh[^\kri3)Yim{-ii0)- (III.ll) 

Im 

The field 4>f3 viewed from the surface E,, is a superposition of solutions to the Helmholtz equation that are regular 
at the origin C„. Using the translation formulas, Eq. (11.30), the field generated by object can be written as 
function of the coordinate x^,, measured from the origin Oa, as 

<^/3(x„) = ikY,Q0,lmY.K^'lmjl'ikra)Yl,m'{^a) • (HI- 12) 

Im I'm' 

This result, in turn, can be substituted into the contribution Spa to the action, leading to the simple result 

S0a[Qa, = I E Qh'm'K^nmQ0,lm + CC. . (III.13) 

Iml'm' 

Note that the contributions to the action that couple fields and sources on different objects make no reference to the 
particular boundary conditions that characterize the Casimir problem. They depend only on the multipole moments 
of the fields and on the geometry through the translation matrix U"''. 



2. Self -interaction terms 
We turn to the terms in 5ci where the field and the source both refer to the same surface, S^: 

SM = \j^ dx(^*(x)(/.,(x) + c.c.) . (III.14) 

For the self-interactions terms, we only use the coordinate system with origin Oa inside the surface E^, and hence 

drop the label a on the coordinates in this section. Since </'q:(x) is continuous across the surface, we can regard the 
in Eq. (III. 14) as the field inside S„, ^in,a) which is a solution to Helmholtz's equation that must be regular at 
the origin Oa, 

a,lmjl 

(fcr)y;„(x). (III.15) 

Im 

Substituting this expansion into Eq. (III. 14), we obtain 

Salea] = \Y1 i'l^»,lmQa,lm + C-C ) , (III. 16) 

Im 

where the Qa.im are the multipole moments of the sources, defined in the previous subsection. 

Finally we relate (pa.im back to the multipole moments of the source to get an action entirely in terms of the Qa,im- 
The field 6 

a,out points outside of 12a obeys Helmholtz's equation and must equal (j)a,in on the surface S. Therefore 
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it can be written as (f>a,in plus a superposition of the regular solutions to the Helmholtz equation that vanish on Sq,, 
as defined in Eq. (11.19), 

0a,out(x) = 0a,in(x) + A<^„(x) = <^„,in(x) + ^ Xa.im {ji{kr)Yi^{yi) + ^ ri?^,i^{k)h!jl\kr)Yi,^, {yi) J . (III.17) 

Im \ I'm' / 

The second term, A(/)q,, vanishes on because T" is the scattering ampHtude for the Dirichlet problem. 

The field we seek is generated in response to the sources and therefore falls exponentially (for k with positive 
imaginary part) as r — > oo. Therefore the terms in Eq. (III. 17) that are proportional to ji{kr) must cancel. Comparing 
Eq. (III. 17) with Eq. (III. 15), we conclude that Xa,im = —^a,iTn, and therefore 

<^a,out(x) = -J2 <t>c.,lm ^ T{^m' lm{k)h^\kr)Yv m' {^) . (III.18) 
Im. I'm' 

On the other hand, 0a,out(x) can be expressed as an integral over the source as in Eq. (III. 2), 

0a,out(x) = ^ (ix'go(x,x',fc)^?a(x'). (III.19) 

Using the partial wave expansion for the free Green's function, Eq. (11.26), we find 

'^a,out(x) =ikY^ Qa,i'm'h'i}\kr)Yi,^,{ii), (III.20) 

I'm' 

and comparing with Eq. (III. 18), we see that 

ikQa,l'm' = ~ %/'m'lm{^)'Pa,lmt 
Im 

or 

<f>a,lm = -ik Yll'^^Vml'm'Qc.^l'm', (III.21) 
I'm' 

where [T«]"i is the inverse of the Dirichlet transition matrix T". When this is combined with Eq. (III. 16), we obtain 
the desired expression for the self-interaction contribution to the action, 

SM = -f E Qhm[T'']Tj,l'm'Qc„l'm' + C.C. . (III.22) 

Iml'm' 



C. Evaluation of the Integral over Sources 

Combining Eq. (III. 22) with Eq. (III. 13), we obtain an expression for the action that is a quadratic functional of the 
multipole moments of the sources on the surfaces. The functional integral Eq. (III. 7) can be evaluated by changing 
variables from the sources, {£>«} to the multipole moments. The functional determinant that results from this change 

of variables can be discarded because it is a common factor which cancels between 5c and Soo- To compute the 
functional integral we analytically continue to imaginary frequency, k — in, k > 0, 

3c(*«) = n / [2^0a^^QyexpJ-^^^Q;[n-io„ + ^^^g;u"''Q/3 + c.c. I , (111.23) 

where we have suppressed the partial wave indices. The functional integral Eq. (III. 23) yields the inverse determinant 
of a matrix M"'' that is composed of the inverse transition matrices [T"]~^ on its diagonal and the translation matrices 
U"/5 on the off-diagonals: 



(III. 24) 
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Finally we substitute into Eq. (11.10) to obtain the Casimir energy, 

^r^i fic r°° , , detMc(m) 

£[C] = — dKln , , { , (III. 25) 

^ ^ n Jo detM^(iK) ^ ' 

where the clctcriniiiaiit is taken with respect to the partial wave iii(Uc:es and the objec;t indic;es a. p. and = 
j-jpaj-ij^^ is the result of removing the objects to infinite separation, where the interaction effects vanish. 
In the special case of two interacting objects Eq. (III. 25) simplifies to 

fir f°° 

£^[C] = — / dKlndet(l - T^V^^T^V^^) , (III.26) 

TI" Jo 

where T", a = 1, 2, and U"'' are the transition and translation matrices for the two objects. 
For three objects £ takes the form, 

£3[C] = — / dK{lndet(l -T^U^iT^Ui^) .TTSuaiTiU^^) 
Jo 

+ Indet [1 - (1 - T3U^^T^U^^)"^(T3U^^ + t3u^1T^U^^)(1 - T^U^1T^U^^)"^(T2U^^ + T^U^^T^U^^)] } . 

(III.27) 

Although this result appears complicated, it admits a simple physical interpretation. The first two terms describe 
scatterings between objects 1 and 2 and objects 1 and 3, and hence correspond to the separate two-body Casimir 
energies of these two pairs of objects. The third term must have a more compHcated form because fluctuation forces 
are not pairwise additive and hence the third term is not simply given by the two-body Casimir energy of objects 2 and 
3. The interaction of objects 2 and 3 involves not only direct scatterings between these two objects but also indirect 
multiple scatterings off object 1. The third term in Eq. (III. 27) contains the resolvent of the operator T^U^-'^T-'^IU"'^'^, 
which can be formally expanded as the series ^^^(T^U'^^T^U^'^)". The resolvent of T^U^^T^U^^ can be treated 
similarly. These series describe multiple scatterings between objects 1 and 3, and objects 1 and 2, respectively, which 
occur as intermediate steps in the scattering processes between objects 2 and 3. Inserting the series into the last 
term of Eq. (III. 27), one can distinguish four qualitatively different scattering processes, described by the following 
operators, 

Oi(m,n) = (T^U^^T^U^^)'"T^U^2 (T^U^^T^U^^)"T^U^^ 
03(m,n) = (WiTiui3-)'"Tr3u3iTriui2(Tr2u2iTiU^2-)"Tr2u23 

04(m,n) = {T^Tn''^W^)"'TVT^Tn^^T^Tn^^T^Tn^yT^\J^^W\ (III.28) 

The operator Oi(0, 0) = T^U^^T^U^^ describes direct scatterings between objects 2 and 3 as if object 1 were absent, 
and hence yields the additive approximation for the energy. All additional operators of Eq. (III. 28) describe the 
non-additivity of the Casimir interaction between three objects. The corresponding scattering processes are depicted 
in Fig. 2. The operator Oi (m,n) describes two scatterings between objects 2 and 3 with m (n) scatterings between 
objects 3 (2) and 1 in between. The operators 02(m, n) {03(171, n)) correspond to a wave that travels from object 3 
(2) to object 2 (3) and is reflected back to object 3 (2) from object 1 after multiple scatterings between objects 1 and 
2, and 1 and 3, respec;tively. Finally, operator 04(171, n) describes only indirect scattering processes between objects 
2 and 3 that all go through object 1. 



IV. APPLICATIONS 



In this section we give a few typical applications of our method. Although the physically interesting case of 
electromagnetism requires vector fields and conducting boundary conditions, the scalar case offers an opportunity to 
check our methods and illustrate their utility. We consider a real scalar field fiuctuating in the space between two 
spheres on which Robin boundary conditions, cp — XadnCp = 0, are imposed. Because a real field has half the oscillation 
modes of a complex field, the Casimir energy in Eq. (III. 26) must be divided by 2, giving 

fir r°° 

£^[C] = — dKlndet(l- W^T^U^i). (IV.l) 

271" Jo 
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2m +1 2(m+l) 



FIG. 2: (Color online) Scattering processes described by the operators of Eq. (III. 28). The directed lines between objects 2 and 
3 describe a wave travelling once between the objects. The zigzag lines correspond to multiple reflections in both directions 
where the number next to the line indicates how many times the wave travels between the objects. To interpret the diagrams 
one starts with a wave at object 3. As an example, we describe explicitly the diagram for the operator Oi (m,n). The wave 
travels from object 3 to object 2, is then reflected 2n times back and forth between objects 2 and 1, travels back from object 
2 to 3, and is finally reflected 2m times back and forth between objects 3 and 1. Each free propagation between the objects is 
described by the translation operator U"'', and each refiection at an object by the transition operator T". 



We allow for different Robin parameters Ai^2 and different radii i?i.2 for the two spheres. This choice allows us to 
study Dirichlet (A/i? — s- 0) and Neumann (A/i? — > oo) boundary conditions on separate spheres as special cases. We 
obtain the Casimir energy as a series in Ri/d and R2/d for large separations d and numerically at all separations. A 
comparison of the two approaches allows us to measure the rate of convergence of our results. We find that for Robin 
boundary conditions the sign of the force depends on the ratios Xa/Ra and on the separation d. We also express 
the Casimir energy between two objects of general shape in a large distance expansion in terms of the coefficients of 
capacitance defined in Sec. TIC. 



A. Interaction of two spheres with Robin boundary conditions: general considerations 

The Robin boundary condition (p — \adn<j) — allows a continuous interpolation between Dirichlet and Neumann 
boundary conditions. Since the radius of the sphere introduces a natural length scale, it is convenient to replace A by 
a dimensionless variable, = Xa/Ra- For > 0, the modulus of the field is suppressed if the surface is approached 
from the outside, while for < it is enhanced. Hence, for negative (a bound surface states can be expected. All 
the information about the shape of the object and the boundary conditions at its surface is provided by the T-matrix. 
For spherically symmetric objects the T-matrix is diagonal and is completely specified by phase shifts Si{k) that do 
not depend on m, 



Iml'' 



,(fc) = 5„,^„„,l(e2^«'«-l) . (IV.2) 



In the discussion of the T-matrix for an individual object we again suppress the label a. The phase shifts for Robin 
boundary conditions are 

where ^ = kR and ji (ni) are spherical Bessel functions of first(second) kind. To apply Eq. (IV. 1), we have to evaluate 
the matrix elements of the transition matrices for imaginary frequencies k — in. Using ji{iz) = i} ^ i: j (2z)1i^yI2(^z) 
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and hi (iz) = —i '-\/2/(7r^)ii';+i/2(-2), we obtain for the T- matrix elements 

T <■ ^ < (i/c+i/2)/i+i/2(^)-^/;+v2(^) ,,,, . 

2 (l/C + 1/2) Ai+l/2(^j - ^^;+i/2(-2) 

where z = kR. 

For two spherical objects we can assume that the center-to-center distance vector is parallel to the z-axis. Then the 
translation matrices simplify to the form presented in Eq. (11.35). For imaginary frequencies the translation matrix 
elements become 

^wi id) = -(-l)"i-''+V(2« + 1)(2Z' + 1) ^ (±1)'"(2/" + 1) I o)i^r'+i/2M), (IV.5) 

where d is the separation distance. 

For a range of (, bound states can appear. This range can be determined by examining the T-matrix at small 
frequencies, where it has the leading form 

where we have suppressed positive numerical coefRcicnts. The amplitudes diverge for ( = X/R = — 1/(Z+1) indicating 
the existence of a pole at k = at this value of (. Indeed, from Eq. (IV. 4) we find that the T-matrix elements have a 
pole in the l^^ partial wave at > if — 1/(Z -|- 1) < C < 0. The values of k; are determined by the equation 

—. = l + l + RKl ' ^)'y (IV.7) 

These poles correspond to bound states for — 1 < C < 0. For any C in this interval, there exists a finite number of 
bound states, which increases as — » 0. In the following, we restrict to ^ > and leave the study of interactions in 
the presence of bound states to a future publication. 

The special case of spheres with Dirichlet boundary conditions has been studied in Ref. [19]. For two spheres 

of equal radius, the matrix ^jp'' A^'™' the notation of Ref. [19] is proportional to our T-'^U^^T^U^^ times 
Ki+i/2{nR)/ Kiij^iI2{k,R)- It is easy to see that this proportionality factor drops out in the final result for the energy 
if one uses Indet = tr In in Eq. (IV. 1) and expands the logarithm around unity. Thus we agree with the results given 
in Ref. [19]. 



B. Asymptotic expansion for large separation 

In this section we consider the Casimir interaction between two spheres due to a scalar field obeying Robin boundary 
conditions, allowing for a different parameter A1.2 on each sphere. The Casimir energy can be developed in an 
asymptotic expansion in R^/d using Indet = tr In in Eq. (IV. 1). Expanding the logarithm in powers of N = 
T^U^^T^U^^, since the T-matrix has no poles in the region of interest we get 




(IV.8) 



We have performed the matrix operations using Mathematica. The scaling of the T-matrix at small k shows that the 
p**^ power of N (corresponding to 2p scatterings) becomes important at order Partial waves of order I start 

to contribute at order if the T-matrix is diagonal in /, which is the case for splieric;ally symmetric; objects. 

Hence the leading terms with p = 1 and / = yield the exact energy to order In the following we will usually 
restrict the expansion to p < 3, Z < 2, yielding the interaction to order d~^. 

1. Equal radii 

We begin with the case Ri = R2 = R- The large distance expansion of the Casimir energy can be written as 




(IV.9) 
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where bj is the coefficient of the term 



d ^ . For general Robin boundary conditions we find 



h = ~ 
hi = - 
h = 



1 



4(1 + Ci)(l + C2)' 
(2 + C1+C2) 



8(1 + Ci)2(l + C2)2' 

- 77 - 231(Ci + C2) - (l88(Cf + Cl) + 625C1C2) - 2 (3(C? + Cl) + 197(CiCl + C2C1)) 
+ 2 (l7(Cf + C2) + 56(C2Cf + CiCl) - 18C?C|) + 2C1C2 (68(C? + Cl) + 155(CiCl + C2C?)) 
+ 2C?Cl (85(C? + Cl) + I24C2C1) + GSdCliCi + C2)l [48(1 + Ci)'(l + 2Ci)(l + (2)^(1 + 20 



(IV. 10) 
(IV.ll) 



6fi = 



- 50 - 175(Ci + C2) - (179(C? + C2) + 592C1C2) - (l8(Ci' + C2) + 565(C?C2 + C2C1)) 

+ (64(Cf + C2) + 6(C?C2 + CiCl) - 460C?Cl) + 2 (ll(Cf + Cl) + 136(Ck2 + C1C2) + 67(C?Cl + CiCD) 
+ 4C1C2 (22(Cf + Cl) + 93(C?C2 + CiCf ) + 52C?Cl) + 2C?Cl (55(C? + Cl) + 92(Cf C2 + CiCD) 



4C?Cl (11(C? + Cl) + 8C1C2) [32(1 + Ci)'(l + 2Ci)(l + C2)'(l + 2C2)]" 



(IV. 12) 



Higher order coefficients can be obtained but arc not shown here in order to save space. For special values of they 
are given below. These coefficients show sonic interesting properties. They all diverge for Ca = —1, where Aq. = —Rai 
and 65 and fog also diverge for Cq = ^1/2- These results are consistent with the emergence of poles in the T-matrix 
at small k when Ca approaches — 1/(^ + 1) and the observation that the partial wave with I = 1 starts to contribute 
only at order . Another interesting property is that some coefficients hj go to zero for Aq, — * 00, which corresponds 
to Neumann boundary conditions. If both Aq- go to infinity, the coefficients bj vanish for j = 1, ... ,6, so that the 
leading term in the Casimir energy is ^ with a negative amplitude. Hence, Neumann boundary conditions lead 
to an attractive Casimir-Polder power law, as is known from electromagnetic field fluctuations. This result can be 
understood from the absence of low-frequency s-waves for Neumann boundary conditions. It is clearly reflected by 
the low frequency expansion of Eq. (IV. 6), which has a vanishing amplitude for Aa ^ ck) if / = 0. If one Aq remains 
finite and the other goes to infinity, only the coefficients 63 and 64 vanish so that the energy scales as with a 
positive amplitude. Since 63 < for Aq > 0, at asymptotic distances the Casimir force is therefore attractive for all 
non-negative finite Aq, and for Aq both infinite. It is repulsive if one Aq is finite and the other infinite, i.e., if one 
sphere obeys Neumann boundary conditions. However, at smaller distances the interaction can change sign depending 
on Aq, as shown below. 

More precisely, one has the following limiting cases. If both Aq = 0, the field obeys Dirichlet conditions at the two 
spheres and the first six coefficients are 



1 

1' 



^4 = -^, &5 



77 



25 



67 



29837 



48 16 2880 

If Neumann conditions are imposed on both surfaces, the coefficients are 



6491 
1152 



63 = 0, 64 = 0, 65 = 0, 66 = 0, 67 = - 



161 , ^ , 3011 , 175 

— , 68 = 0, h, = -— b,, = -—, 



(IV. 13) 



(IV. 14) 



clearly showing that the asymptotic interaction has a Casimir-Polder power law ~ 0{d~'^). Also, as in the electro- 
magnetic case, the next to leading order 0{d~^) vanishes [1]. Therefore we have included the two next terms of the 
series. If A2 ^ 00 (Neumann conditions) and Ai remains finite, we have 



63 = 0, 64 = 0, 65 = 



17 



11 



48(1 + Ci)' 



32(1 + Ci)2 



67 = 

68 = 



1989 + 4736Ci + 1895C? - 2192C? - I6IOC1 

480(1 + Ci)3(l + 2Ci) 
5(94 + 143Ci + 76Cf - 45Cf ) 



288(1 + Ci)^ 

Thus, for mixed Dirichlet /Neumann boundary conditions, the result is 

17 11 

63 = 0, 64 = 0, 65 = — , 66 = — , 67 = 



663 
160' 



235 
144 



(IV. 15) 



(IV. 16) 
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Finally, if A2 = (Dirichlet conditions) and Ai is finite, we obtain 



4(1 + Ci)' " 8(1 + Ci)2' 
77 + 23IC1 + I88C1 + 6C1 - 34C^ 



65 = - 

ifl -I- 1 7f,r, 4- 1 7Q/-? 4- 1 8^3 _ fi4^4 _ 99^5 

(IV. 17) 



48(1 + Ci)3(l + 2Ci) 
50 + 175Ci + 179Ci + 18C? - 64Cf - 22Cf 
32(l + Ci)4(l + 2Ci) 



where we have not shown higher order coefficients. It is important to note that the series in Eq. (IV.9) is an asymptotic 
series and therefore cannot be used to obtain the interaction at short distances. 



2. Unequal radii 

To study the individual contributions from two objects to the various terms in the large distance expansion, it is 
instructive to study two spheres of different radii i?i and R2. For simplicity, we focus on Dirichlet and Neumann 
boundary conditions. By expanding the energy in powers of N as before, we get an asymptotic expansion of the 
Casimir energy, 

^-?i£M.)(f)"". (IV.18) 

j=3 ^ ^ 

where the coefficients now also depend on r] = R2/R1. For Dirichlet boundary conditions on both spheres, the 
coefficients are 

- r] ~ ri + if ~ 34(77 + r?3) + 9^2 2(?? + r;^) + 23(7?^ + r/^) 

63 = -4, 64 = — 65 = , h = ' ^^^-^^^ 

- _ _ 8352(?7 + 77^) + 1995(772 + 77*) + 3898077^ h - - "^^^^(^ + '^^) + 5478(77^ + 7;^) + 2357(77=^ + 7?^) 

^ ~ 5760 ' ^~ 2304 ' 

while for Neumann boundary conditions on both spheres, we have 

-^.=-"1^, ~^8=o, 6,=-«^, w,=- ''%:^'\ (IV.20) 

96 384 256 ^ ' 

and for Dirichlet conditions on the sphere with radius R\ and Neumann conditions on the one with radius ii2, we 
obtain 

- 17773 - II773 - 77^(1610 + 37977^) - 577^(24 + 6777^ + 3773) 

^^ = ^' ^'=^^ 480 ' 288 • (^^-''^ 

For like boundary conditions, the coefficients satisfy the symmetry condition r]^~^bj{l/r]) = bj{rj). The coefficient 
63 for Dirichlet conditions becomes positive if 77 is sufficiently large or small compared to unity. Due to the absence 
of monopoles, Neumann conditions on one or both spheres lead to coefficients that scale at least as rf . For the 
same reason, not all powers of 77 contribute to the coefficients if Neumann conditions axe imposed. Note that the 
interaction between a sphere and a plate cannot be obtained from this asymptotic expansion as the limit 77 — > 
since the expansion assTimes that Ri, R2 ^ d. However, the sphere-plate interaction can be computed by the same 
technique employed here if it is applied to the Green's function of the semi-infinite space bounded by the plate instead 
of the free Green's function. 



C. Numerical results for Robin boundary conditions on two spheres at all separations 

The primary application of our analysis is to compute the Casimir energy and force to high accuracy over a broad 
range of distances. However, to obtain the interaction at all distances, Eq. (IV. 1) has to be evaluated numerically. We 
shall see that the domain where our method is least accurate is when the two surfaces approach one another. That is 
the regime where semiclassical methods like the proximity force approximation (PFA) become exact. Because of its 
role in this limit, and because it is often used (with little justification) over wide ranges of separations, it is important 
to compare our calculations with the PFA predictions. 
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1. PFA for Robin boundary conditions 



In the proximity force approximation, the energy is obtained as an integral over infinitesimal parallel surface 

elements at their local distance L, measured perpendicular to a surface S that can be one of the two surfaces of the 
objects, or an auxiliary surface placed between the objects. The PFA approximation for the energy is then given by 



^PFA = £\\{L)dS , 



(IV. 22) 



where £\\{L)/A is the energy per area for two parallel plates with distance L. Functional integral techniques can be 
used to compute the energy for parallel plates with mixed (Robin) boundary conditions [10]. For brevity we only 
quote the results [29]. The energy is 



£:||(L) = |^$(Ai/L,A2/i), 
and the force per area is 



he 



$(Ai/i,A2/L) 



du In 



1 



(l-uAi/L)(l 



UA2/L)^_2„ 



^11 (L) = j^^'{X,/L,X2/L), ^'{X^/L,X2/L) 



dui 



(l + uAi/L)(l + uA2/L) 



{l + uXi/L){l + uX2/L) 
{l-uXi/L){l-uX2/L)' 



(IV. 23) 



(IV. 24) 



The amplitudes ^{Xi/L, X2/L), ^'{Xi/L, X2/L) are finite for all Xa/L > 0. In general, the integrals in Eqs. (IV.23), 
(IV. 24) have to be computed numerically. The PFA predictions as functions of Ai,2 have considerable structure, which 
we illustrate in Fig. 3. 

The behavior of the PFA at asymptotically small or large A„/L determines the Casimir interaction as L ^ 0. For 
all non-zero values of Ai,2, we take Xa/L 00, but for the Dirichlet case, A = 0, the limit Xa/L applies. For 
parallel plates with Robin boundary conditions, in the limit Xi^2/L » 1 we obtain the result for Neumann boundary 
conditions on both plates. 



$(Ai/L,A2/L)^$c 



1440 



(IV. 25) 



and for X\^2/L <C 1 we obtain the identical result for plates with Dirichlet conditions. Finally for X\^2/L <C 1 <C X2,i/L, 
we obtain the parallel plate result for unlike (Dirichlet/Neumann) boundary conditions. 



$(Ai/X,A2/L) 



*o = 



11520 



(IV. 26) 



The last case is relevant at short distances if one of the A^ = 0. For two spheres of radius R and center-to-center 
separation d with Robin boundary conditions, the PFA results can now be obtained easily from Eq. (IV. 22). In terms 
of the surface-to-surface distance L = d — 2R, we get 



±. 77 Rhc 

° 2{d- 2RY ' 



(IV. 27) 



where the -I- applies if one and only one Aq, = 0, and the — in all other cases. Hence, at small separation the interaction 
becomes independent of A^, in the sense that it only depends on whether one Aq, is zero. In fact, the form of the 
T-matrix of Eq. (IV. 4) shows that the energy becomes independent of A^ for d ^ min(Ai, A2). Indeed, if we write the 
frequency as k = u/d, the translation matrix depends only on u, and the T-matrix can be written as 



Imlri 



(-1) 



^ (d/Xa + d/2R^)Ii+^/2{uRa/d) 



■"-^i+l/2("^"/'^) 

2 (d/A„ + d/2R„)Ki+i/2{uRa/d) - uK[^^^^{uRa/ d) 



(IV. 28) 



so that the energy depends only on d/Ra for d <C A^. This universal behavior is confirmed by a numerical evaluation 
of the energy that is described below. 

With the results obtained above, we can analyze the sign of the interaction between plates and spheres at both 
asymptotically large and small distances. Since the PFA result is expected to hold in the limit where the distance 
tends to zero, Eq. (IV. 27) predicts the sign of the interaction between spheres in the limit of vanishing distance. 
In the limit of large distances, we can compare the results for parallel planes from Eqs. (IV. 25) and (IV. 26) to our 
calculations for two spheres. We find that the sign of the asymptotic interaction depends on the choice for Xa and 
is identical for plates and spheres. Hence, we obtain a complete characterization of the sign of the interaction at 
asymptotically large and small distances for the plate and sphere geometry, which is summarized in Table I. However, 
as we have seen above, the power law decay at large distance is quite different for plates and spheres. 
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(for plates: A1/A2 or A2/A1 > 2.8) 
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]0,oo[ 


00 




+ 





TABLE I: The sign of the Casimir force between two plates and two spheres with Robin boundary conditions at asymptotically 
small and large surface-to-surface distance L. The sign in these two limits is identical for plates and spheres. Here "— " and 
indicate attractive and repulsive forces, respectively. 




1-hAi/L 

FIG. 3: (Color onhne) Amplitudes $(Ai/L, A2/L) of the Casimir energy (from eq. (IV.23)) and $'(Ai/L, A2/I/) of the Casimir 
force (from eq. (IV. 24)) for two parallel plates of distance L with Robin boundary conditions. L = Q corresponds to the 
point (1, 1) unless one or both plates have Dirichlet (A^ = 0) boundary conditions. The color coding corresponds to the 
magnitude of the energy amplitude with blue corresponding to positive energy and red to negative. The dashed (black) curves 
correspond to zero energy while the solid (red) curves represent zero force. The curve running from the origin to (1, 1) shows 
as an example the interaction energy for A2/A1 = 5 as function of the distance L. For finite A^ the interaction is always 
attractive at asymptotically large and small distances but can be repulsive at intermediate distances if the ratio A1/A2 > 2.8 
(or A2/A1 > 2.8). For Ai = A2 the interaction is always attractive. The amplitudes are independent of L for the special cases 
where the Aq, are either zero or infinite: The interaction is attractive (negative energy) for Ai = A2 = 0, Ai = A2 = cxd and 
repulsive (positive energy) for Ai = 0, A2 = 00 and Ai = 00, A2 = 0. 



2. Casimir forces for all separations 

To go beyond the analytic large distance expansion, we compute numerically the interaction between two spheres 
of the same radius R with Robin boundary conditions. Guided by the classification of the Casimir force according to 
its sign at small, intermediate and large separations, we discuss the six different cases listed in Table I. Our numerical 
approach starts from Eq. (IV. 1). Using the matrix elements of Eq. (IV. 4) and Eq. (IV. 5), we compute the determinant 
and the integral over imaginary frequency numerically. We truncate the matrices at a finite multipole order I so that 
they have dimension (1 -I- Z)^ x (1 -I- Z)^, yielding a series of estimates for the Casimir energy. 

£(1) gives the exact result for asymptotically large separations, while for decreasing separations an increasing 
number of multipoles has to be included. The exact Casimir energy at all separations is obtained by extrapolating 
the series {£^^^} to I — > 00. We observe an exponentially fast convergence as IS a 

constant of order unity. Hence, as the surfaces approach each other for d 2R, the rate of convergence tends to 
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zero. However, we find that the first I = 20 elements of the series are sufficient to obtain accurate results for the 
energy at a separation with R/d — 0.48, corresponding to a surface-to-surface distance of the spheres of L = 0.083i?, 
i.e., approximately 4% of the sphere diameter. In principle our approach can be extended to even smaller separations 
by including higher order multipoles. However, at such small separations semiclassical approximations like the PFA 
start to become accurate and can be also used. 

The results for Dirichlet and Neumann boundary conditions are shown in Fig. 4. All energies are divided by EppA, 
given in Eq. (IV. 27), with the corresponding amplitude $o (repulsive at small separations) or $q (attractive at small 
separations). For like boundary conditions, either Dirichlet or Neumann, the interaction is attractive at all separations, 
but for unlike boundary conditions it is repulsive. At large separations the numerical results show excellent agreement 
with the asymptotic expansion derived above. Note that the reduction of the energy compared to the PFA estimate 
at large distances depends strongly on the boundary conditions, showing the different power laws at asymptotically 
large separations. In the limit of a vanishing surface-to-surface distance {R/d 1/2), the energy approaches the PFA 
estimate in all cases. Generically, the PFA overestimates the energy: fppA is approached from below for R/d — > 1/2, 
except in the case of Dirichlet boundary conditions on both spheres, where the PFA underestimates the actual energy 
in a range of 0.3 < R/d < 1/2. The deviations from the PFA are most pronounced for Neumann boundary conditions. 
At a surface-to-surface distance of L = 3R {R/d = 0.2), the PFA overestimates the energy by a factor of 100. 

Casimir interactions for Robin boundary conditions with finite are shown in Fig. 5. If Ai = A2 the interaction 
is always attractive. If the A^ are not equal and their ratio is sufficiently large, the Casimir force changes sign either 
once or twice. This behavior resembles the interaction of two plates with Robin boundary conditions. However, the 
criterion for the existence of sign changes in the force now depends not only on A1/A2, but on both quantities Xi/R 
and X2/R separately. Even with A1/A2 fixed, for smaller Xa/R there can be sign changes in the force, while for 
larger A„/i? the force is attractive at all distances. When the ratio A1/A2 is sufficiently large (or formally infinite for 
Dirichlet or Neumann boundary conditions), we can identify three different generic cases where sign changes in the 
force occur: 

• First, we consider Dirichlet boundary conditions (Ai = 0) on one sphere and a finite non- vanishing X2/R at 
the other sphere. Figure 5(a) displays the energy for A2/i? = 10 as a typical example. At large distances the 
energy is negative, while it is positive at short separations with one sign change in between. The asymptotic 
expansion of Eq. (IV. 9) with the coefficients of Eq. (IV. 17) yields the exact energy at separations well below the 
sign change. While the expansion predicts qualitatively the correct overall behavior of the energy, it does not 
yield the actual position of the sign change correctly. Of course, for the Casimir interaction between compact 
objects, the sign of the force = —d£/dd is the physically important quantity, not the energy. The distance at 
which the force vanishes cannot be deduced directly from the slope of the curve for £ = E/S-p-pA, since one has 



£'{d) ^ 



£pFx{d) 

The force vanishes at the distance do if £'{dQ) = 0, so that 



(IV.29) 



^'(rfo) = -^[^^^ido) . (IV.30) 

Hence the distance at which the force vanishes is determined by the position do where the curve of the auxiliary 
function t{d) = T{d/R — 2)^ is tangent to the curve of £. The two unknown quantities do and r are then 
determined by the conditions £{do) = t{do) and £'{do) = t'(do)- This procedure allows us to obtain the distance 
at which the force vanishes easily, without computing derivatives numerically. The tangent segment of the curve 
for t{d) is shown in Fig. 5(a) as a dotted line. From this construction we find that at a distance the force 

changes from attractive to repulsive for decreasing separations. The position corresponds to a niininuim 

of the energy and decreases with decreasing X2/R, so that in the limit X2/R^ it approaches the case of two 
spheres with Dirichlet boundary conditions, where the force is always attractive. 

Second, we study Neumann boundary conditions on one sphere and a inite non- vanishing A2/i? at the other 
sphere. As an example we choose again X2/R = 10, as shown in Fig. 5(b). The energy is positive at large 
distances and becomes negative at small distances. The asymptotic expansion with the coefficients of Eq. (IV. 15) 
is found to be valid well below the separation where the sign of the energy changes. Hence, the expansion 
describes the behavior of the energy qualitatively, but does not predict the precise position of the sign change. 
The sign change of the force can be obtained by the method described above. At a position d+=^_, the force 
changes from repulsive to attractive with decreasing separation and the energy is maximal. A decreasing 
(increasing) A2/i? shifts to smaller (larger) separations. This result is consistent with an entirely repulsive 

(attractive) force for Neumann-Dirichlet (Neumann-Neumann) boundary conditions. 
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FIG. 4: (Color online) Casimir energy for two spheres of radius R and center-to-center distance d: (a) Dirichlet boundary 
conditions for both spheres, (b) Neumann boundary conditions for both spheres, (c) Spheres with different boundary conditions 
(one Dirichlet, one Neumann). The energy is scaled by the PFA estimate of Eq. (IV. 27). The solid curves are obtained by 
extrapolation to / ^ oo. For the smallest separation, the extrapolation uncertainty is maximal and indicated by an error bar. 
The dashed curves represent the asymptotic large distance expansion given in Eq. (IV. 9) with the coefficients of Eqs. (IV. 13), 
(IV. 14) and (IV. 16), respectively. 

• The third case is obtained if both Aq, are finite and non-zero. A typical example with Xi/R — 20 and \2/R = 1 is 
shown in Fig. 5(c). The energy is negative both at large and small separations but turns positive at intermediate 
distances. The asymptotic expansion applies again at sufficiently large separations beyond the position where 
the energy becomes positive. For values of the ratio A1/A2 that are larger than an i?-dependent threshold, 
the force changes sign twice, so that it is repulsive between the separations and The energy has 

a minimum (maximum) at d-=>+ (d+=^-). If Ai/i? increases and X2/R decreases, the repulsive region grows 
until eventually the force becomes repulsive at all separations, corresponding to the limit of Dirichlet/Neumann 
boundary conditions. Decreasing Xi/R and increasing X2/R reduces the interval with repulsion. In this case, 
first the zeros of the energy disappear, leaving negative energy at all distances but still a repulsive region, and 
then the two positions where the force vanishes merge, leaving an entirely attractive force. 
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FIG. 5: (Color online) The Casimir energy for two spheres with different Robin boundary conditions for finite Aq.: (a) Dirichlet 
boundary conditions and \/ R = 10, (b) Neumann boundary conditions and A/71 = 10, (c) Xi/R = 10 and A2/-R = 1. The solid 
curves correspond to extrapolated results for / ^ 00, and the dashed curves represent the asymptotic large distance expansion 
given in Eq. (IV. 9) with the coefficients of Eqs. (IV. 17), (IV. 15), (IV. 10), respectively. For logarithmic plotting, the modulus 
of the energy is shown, and the sign of the energy is indicated at the bottom. The range of separations with a repulsive force is 
shaded. The points of vanishing force occur where an auxiliary function (dotted curves) is tangent to the solid curve, see text 
for details. 



D. Interaction in terms of low frequency scattering data 



We noted in Sec. II C that in the limit of vanishing frequency, the S-matrix is related to tensor generalizations 
of the capacitance. At higher frequencies, this tensor receives corrections that can be expanded in frequency. The 
object's shape is thus encoded in the tensor expansion coefficients of the S-matrix. We can use these results to obtain 
the Casimir interaction between non-spherical objects in a large distance expansion. From such an analysis one can 
obtain information about which pieces of the S-matrix contribute to the Casimir interaction at a given order in inverse 
separation. 

In general, according to Eq. (11.29), the T-matrix of object a can be written in a small K-expansion of the form 

■i+i' °° 

■Tm'lm (2f I 1)!!(2^' 1)!! , Q\l' m' Im^ ■ i^-l v .u-l; 
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Due to the symmetry of the T-matrix, the coefficients obey the relation 



- ^ 

q\L'm'Lm 



^ q;lml' m' 



(IV.32) 



We assume again that the origins of the two objects are located on the 2:-axis at a separation d. Then we can use 
the translation matrix elements of Eq. (IV. 5). The Casimir energy can be obtained from the 1ow-k expression for 
the T-matrix of Eq. (IV.31) by a simultaneous expansion in the number of partial waves, the number of scatterings 
(powers of N), and the imaginary frequency k. The analysis based on the general low- frequency structure of the 
T-matrix for objects of arbitrary shape and boundary condition becomes increasingly complicated with decreasing 
separation between the objects. Since the T-matrix cannot be assumed to be diagonal in /. partial waves of order I 
can contribute to the energy at order (for T-matrices that are diagonal in I, they can only contribute at order 

^-(3+2;)^ We consider here only partial waves up to order I = 3 and four scatterings, which allows us to obtain the 
energy to order d~^. The energy can again be written in the form 



he Bj 
J=3 



(IV.33) 



The coefficients Bj are completely specified by the low-frequency behavior of the T-matrix, i.e., they can be obtained 

from the coefficients C^.i,^,i^ in Eq. (IV.31). In the following, we assume that the origin of the object's coordinate 
frame has been chosen such that the dipole response to an applied constant potential vanishes, i.e., C^i^oo — 
Cg.ooim = for m = —1,0, 1. Using the symmetry relation of Eq. (IV.32), the result for the coefiicients can be written 
as 



Bi 
B5 



'^0;0000'^0;0000> 



('^0;0000^1;0000 + ^0;0000^1;i 



0000 y > 



1 

16 

5 



[Co;0000^0;OOOo]^ 



0000^0:1010 



QiiOoooQiaoio) 



(C'o;002oC'o;0000 " 



(IV.34) 

(IV.35) 
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'32 
7 /5 



(C'o;0000^0;1020 



^0;0000^0;102o) 



C'o;OOOoC'2;0000 



■ C'o;0000^2;OOOo) ' (IV. 36) 
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'l6 
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(>'-^0;0000'-^2;0010 ~ '-^0;0000'-^2;0010 '-^1;0000'-^1;0010 ~ '-^l;0000'-^l;0010j 

+ ^1; 0000 ^2; 0000 + ^0;0000^3;0000 + C'o:OOOoC'3;OOOo) • 



(C'l:0000^2;0000 



(IV.37) 



This result deserves several comments. First, the above expressions for the coefficients have been substantially 
simplified by the requirement that Cf^oo = ^ooim = 0) which eliminates 2, 5 and 22 terms from B4, B5 and B^^ 
respectively. Most contributions originate from two scatterings, which can be seen from the appearance of a product 
of two coefiicients C'^j^i'm'- O^ily the leading terms of B^ and B^ are composed of a product of four coefiicients 
^q;imi'm' ^^'^ heucc rcsult from four scatterings. The latter terms are completely determined by the capacitance of 
the objects and its finite frequency corrections, C^oooo- "^^^^ result clearly shows the relevant number of partial waves, 
which increases by one with each additional power of l/d. At leading order, ~ l/d'^, the coefiicient B3 is given by 
the product of the capacitances Co^oooo objects. At the next order, contributions from I = 1 partial waves are 

absent due to the choice of origin. At higher orders, the coefficients Bj contain contributions from partial waves up 
to ^ = j — 3. Note also that not all terms are symmetric under an exchange of the two objects, C^_;^;,^, ^ C^q^imi'm" 
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The terms for which I + 1' is an odd number change sign under an exchange of the objects. This is related to the fact 
that the expansion of the T-matrix in Eq. (IV. 31) assumes that the local coordinate systems of the two objects have 
the same orientation. Hence, the objects "see" each other along different directions, i.e., along the direction of the 
positive and negative 2-axis, respectively. Up to coefficient Bq, all terms have m = m' = 0, which is again due to the 
choice of origins eliminating all terms with nonzero m up to the order considered here. 

The above result can be applied to two objects of arbitrary shape and orientation, since the direction connecting 
the origins of the objects can be always defined as the z-axis. The orientation of the objects with respect to this 
axis is then specified by the coefficients C^;^;/^/, which depend on the angle that a reference direction, fixed to the 
object, forms with the z-axis. Hence, the dependence of the Casimir interaction on the orientation of the objects can 
be obtained from this general result by computing the low-frequency expansion of the T-matrix elements for objects 
with an arbitrary orientation with respect to the 2;-axis. Since the capacitance coefficients C^oooo independent 
of orientation, a dependence on orientation can occur only at order d'^^. A more detailed study of the orientation 
dependence and the reduction of the number of independent expansion coefficients due to symmetries of the objects 
is left to a future publication. 
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Throughout the paper we considered Dirichlet boundary conditions, ^ = 0. In fact, our final result, Eq. (III. 25), 
which expresses the interaction Casimir energy of many compact objects in terms of the transition and translation 
matrices, is independent of the choice of boundary conditions, provided the T-matrix is the one appropriate to the 
boundary conditions of interest. Here we show this result for the Neumann case. The general Robin case is left to 
the reader. 

Neumann boundary conditions are implemented by replacing 4>(x) by 9„(/)(x) in Eq. (11.12), leading to an expression 
for the partition function analogous to Eq. (11.13) with 0(x) 9„0(x). Like the Dirichlet case, this case also has an 
analogy to electrostatics, namely a complex field coupled to a set of surface dipole densities ga{x). By analogy, in 
this case (?„(^ will be continuous throughout space, but itself will jump by Qa{^) across the surface S^- Therefore 
the classical equations of motion analogous to Eq. (IH.l) are 



and the normal derivative of the free Green's function generates the field associated with these dipole sources, 



V. 
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APPENDIX A: NEUMANN BOUNDARY CONDITIONS 



(V' + A;2) 0ei(x) = 0, forx^S^, 

A(/)ci(x) = -£'a(x), for X e E, 
A9„<?ici|x = o> for X e , 



(A.1) 




(A.2) 



where 9„' denotes the normal derivative at the point x' and acts only on Qq. 
The evaluation of the classical action proceeds in analogy to the Dirichlet case. 




(A.3) 



First we evaluate the terms with a ^ /3, 




(A.4) 
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In analogy to Eq. (III.9), 

which serves to define multipole moments of the dipole layer density, 

P0,im = / dx0dn [ji{kr0)Yil^{x0)] Qoixf)) . (A.6) 

As in the Dirichlet case, the Hankel functions defined with respect to Op can be expanded in terms of Bessel functions 
about Oa with the help of translation formulas. Substituting into Eq. (A. 4) leads to a result analogous to Eq. (III. 13), 

.S/3a[Pa, = ^ ^ Ph'^'^i'lmPp^l^ + ^-C • (A.7) 
Iml'm' 

When a = /?; wc also proceed in analogy to the Dirichlet case. The interior field is defined as in Eq. (III. 15) and 
substitution into Eq. (A. 3) leads to 

[Qa] = \Y^ {<t><,,lmPa,lm + C-C ) , (A.8) 

Im 

in analogy to Eq. (III. 16). As in the Dirichlet case, we determine the field expansion coefiicients, 4>a,imi by constructing 
the exterior field (t>out,a in two different ways. First we construct it following the logic of Eqs. (III.17)-(III.20), but 

now with the aid of the Neumann T-matrix , since the normal derivative of <f>a must be continuous across Eq, . Then 
we construct (t>ont,a directly from the Neumann multipole moments using Eq. (A. 2). Equating the two we obtain 

Km = -ik Y.i'^Trlvm'Pc.^l'm', (A.9) 
I'm' 

where T" is the Neumann T-matrix. Substituting into Eq. (A.8) we obtain the final result, 

Sa[Pa] = -y E Phm[TTm,m'Pc^,l'm' + CC. , (A.IO) 
Ijn 

analogous to Eq. (III. 22). The rest of the derivation follows the Dirichlet case closely. 



APPENDIX B: CASIMIR ENERGY AND FUNCTIONAL INTEGRALS 



Here we show the equivalence between our functional integral formalism and the standard representation of the 
Casimir energy in terms of zero-point energies. We consider a real scalar field, (p{x, t), that obeys the wave equation 
in a domain V, and obeys a Dirichlet boundary condition on E, the boundary of V, which need not be connected. 
Let 0a(x) and —e'l^/h'^c? be the eigenfunctions and eigenvalues of the Laplacian in D, 

-V2</.„(x) = ^0a(x) in V, 

K{x) =0 on E. (B.l) 

Let — „/fi^c^ be the corresponding eigenvalues when the objects arc moved to arbitrarily large separation from 
one another. For simplicity we assume that the spectrum is discrete in both cases. 

Consider the Minkowski space functional integral over the field ip{x,t), as introduced in Eq. (II. 1) but now for a 
real field, 

Z[C] = J [I?v']cexp|^^rfi^dx(^l(9t^)2-(V^)2^| . (B.2) 

By definition, only periodic paths, (p{x,0) = (p{'x.,T), that obey the Dirichlet boundary conditions axe included. Since 
(p(x, t) is periodic in T, it can be expanded in a Fourier series, 

oo 

<p{x,t)= E <Pn(x)e-2™*/T^ 
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with <y9_„(x) = 93* (x). Because the boundary conditions are time independent, the functional integral over each 
frequency mode of can be done independently. Therefore Z[C\ can be written as an infinite product, 

Z[C]= n /[2?¥'.]cexp|f ^dx(^(^^)Vnr-|V<^n|'j| . (B.3) 

At the end of this calculation we will divide by Z^ooj the functional integral defined when all the objects are removed to 
arbitrarily large separation. Therefore we ignore all multiplicative factors in Z\C\ that are independent of the location 
of the constraining surfaces, beginning with the Jacobian generated by the change of variables from (^(x, f) to c/5„(x). 

The real eigenfunctions, {0a (x)}, defined in Eq. (B.l) form a complete orthonormal set of functions on V, so the 
field i^n(x) can be expanded as 

a 

where {c„q} are complex numbers with C-na = c*^ . When this expansion is substituted into Z[C], the functional 
integral reduces to an infinite product of ordinary integrals over the coefficients {c„q,}j 

m = ft n/ E f^- (^)l M'] ' (B-4) 

n=—oo a \ ot \ / ) 

where we have dropped factors common to Z[C] and Zrxi- The integrals over the {cna\ can be performed (taking T 
complex with a negative imaginary part, and analytically continuing back to T real at the end), leading to 



-1/2 

(B.5) 



n=—oo a 



Then, interchanging the products over n and a, Z[C] can be written, 

(B.6) 



£2 
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where we have separated out the n = term and combined each positive and negative n into a single term. Factor- 
ing out the configuration independent factor (^^) ^^'^ using the infinite product representation of the sine 
function. 



sma; 



n 

n=l 



we obtain 

z[c] = n 



e 2R 

— e 



(B.7) 



Finally we define T = —iA/c, send A — > 00, divide by Z^o, and take the logarithm, obtaining 



lim In iZ[C]/Z^) = -^y^iea - ea.oo) , (B.8) 
A— >c!o 2nc '-^ 

which, using Eq. (II. 4), gives the standard expression for the Casimir energy of a real scalar field, 

a 

Note that if we had set T = and not taken the T — > c» limit, we would have obtained the partition function for 

a Bose "gas" of scalar particles occupying the eigenstates with energy £q,) the starting point for a quantum statistical 
mechanical treatment of the Casimir effect at temperature 1 //3. 
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APPENDIX C: TRANSLATION FORMULAS 

For completeness, we sketch the derivations of the translation formulas quoted in Eqs. (11.30) and (11.31) with /? = 1 
and a = 2. See Fig. 1 for the configuration of the two objects. 

1- "l^/'m'im (X12) 

First consider the equation relating the regular partial wave solution about X2 to those expanded about O^. Since 
X2 = X21 + xi we can write 



gikx2 _ gikX2igikxi 



for an arbitrary vector k. If we use the familiar partial wave expansion, 

e^kx ^ 47r^z'j;(fcr)y,;;,(k)ll„(x), 



(C.l) 



(C.2) 



in all three terms, multiply by the factor y/m(k) and integrate over k, we obtain 
ji(fc^2)l;„(x2) = 47r ^ ( j dky,„(k)y;„,(k)y;?,„„(k)] 3i"{kd^2)3v{kn)Yv'm"{±2i)Yvm'{^i) ■ (C.3) 

l'm'l"m" ^ ^ 

The integral over three spherical harmonics is [30] 

I dky,„(k)r;„,(k)r;,„„(k) = (-i)™'^'^^^^^. r r ) j^, _r ^ (c.4) 

where X = 2l + 1. Substituting into Eq. (C.3) and regrouping terms we find. 



ji{kr2)Yim{yi2) = XI 



\/47rAA^/-'(-l)'"'jHfcri)y;w(xi) • 



(C.5) 



Comparing this expression with the definition of the translation matrix, V, given in Eq. (11.30), we can read off the 
matrix elements V/^;,„,(Xi2), 

V/4';™(Xi2) = ^/^i^'-'i-ir' E *'"(-l)""v^ '0 'o){L -i' Jl")Mkdi2)Y,,m"{^2i) • (C.6) 

l"m" 

Finally we use X21 = — X12 and Fp/m"(— X) = (—1)' Yi"„i»{X.), together with the fact the the 3j-symbol restricts 
I' + I" — I to be even and m' + m" to equal m, all of which yields 

V/4'^„^(Xl2) = s/^i'-'\-ir E '0 ) (™ -i' Z")jl"{kd^2)Yl.m"{±12) , (C.7) 

l"m" 

the analogue of Eq. (11.31). 

We begin by equating the partial wave expansion (with respect to O2) of the free Green's function to its plane wave 
representation. 



e^o(x,X2,fc) = ik^ji{kr)hf'{kr2)Yi*^{x)Yi^{i2) 



dq e*i(''=^-'') 
(27r)3 g2_fc2 • 



(C.E 
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Without loss of generality we take |x2| > |x| with respect to 02- We project out the (/m)*^ partial wave by multiplying 
by Yimi^i-) and integrating over x, 



ikji{kr)hl^\kr2)Yim{i2) = J dmmi^) J 



(27r)3 q^-k^ ■ 



(C.9) 



Next we write e*'*^''^^''-' = e*i(-^2i+xi-x) ^^^^ expand each exponential in partial waves using Eq. (C.2). The x integral 
can be performed using the orthonormality of the Yim', the q integral is of the form of Eq. (C.4). The substitution of 
X21 = — X12 proceeds as in the previous case and we obtain 



ikji{kr)h'l^\kr2)Yimiyi2) 



i-iy 



47r 



/ J \0 / \m —m —m j 



'm'V'm" 



f 

.Jo 



dq 



(p' — k"^ 



This integral can be found in Ref. [31], 



dq- 



ji"{qdi2)ji'{qri)ji{qr) 



jji"{qdi2)ji'{qri)ji{qr) = '^ji{kr)ji'{kri)h[},\kdi2). 



Yi„m"{^12)Yi,m'{^l). 



(C.IO) 



Substituting this result into Eq. (C.IO), the form of ^^;,^,(Xi2) quoted in Eq. (11.31) follows straightforwardly. 
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